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Abstract 

We reexamine remarkable connection, first discovered by Beuk- 
ers, Kolk and Calabi, between C(2n), the value of the Riemann zeta- 
function at even positive integer, and the volume of some 2n-dimensio- 
nal polytope. It can be shown that this volume equals to the trace 
of some compact self-adjoint operator. We provide an explicit ex- 
pression for the kernel of this operator in terms of Euler polynomials. 
This explicit expression makes it easy to calculate the volume of the 
polytope and hence (,{2n). In the case of odd positive integers, the 
expression for the kernel enables to rediscover an integral representa- 
tion for Q(2n + 1), obtained originally by different method by Cvijovic 
and Klinowski. Finally, we indicate that the origin of the Beukers- 
Kolk-Calabi's miraculous change of variables in the multidimensional 
integral, which is at the heart of all of this business, can be traced 
down to the amoeba associated with the certain Laurent polynomial. 

1 Introduction 

In a nice little book [1] Vladimir Arnold has collected 77 mathematical prob- 
lems for Kids from 5 to 15 to stimulate the development of a culture of 
critical thinking in pupils. Problem 51 in this book asks to calculate the sum 
of inverse squares of the positive integers and prove the Euler's celebrated 
formula 
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Well, there are many ways to do this (see, for example, [21 El lU El El E] 
and references therein), some maybe even accessible for kids under fifteen. 
However, in this note we concentrate on the approach of Beukers, Kolk and 
Calabi (8], further elaborated by Elkies in [9j. This approach incorporates 
pleasant features which all the kids (and even some adults) adore: simplicity, 
magic and the depth that allows to go beyond the particular case ([I]). The 
simplicity, however, is not everywhere explicit in [Sj and [9], while the magic 
longs for explanation after the first admiration fades away. Below we will try 
to enhance the simplicity of the approach and somewhat uncover the secret 
of magic. 

The paper is organized as follows. In the first two sections we reconsider 
the evaluation of £(2) and £(3) in order technical details of the general case 
not to obscure the simple underlying ideas. Then we elaborate the general 
case and provide the main result of this work, the formula for the kernel 
which allows to simplify considerably the evaluation of ({2n) from [HI |9j 
and re-derive Cvijovic and Klinowski's integral representation [10] for ((2n + 
1). Finally, we ponder over the mysterious relations between the sums of 
generalized harmonic series and amoebas, first indicated by Passare in [7]. 
This relation enables to somewhat uncover the origin of the Beukers-Kolk- 
Calabi's highly non-trivial change of variables. 



The sum (fT]) is just ((2) which we will now evaluate following the method 
of Beukers, Kolk and Calabi [8] . Our starting point will be the dilogarithm 
function 



2 Evaluation of ((2) 



Recall the definition of the Riemann zeta function 
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Clearly, Li 2 (0) = and Li 2 (l) = C(2). Differentiating (j3J), we get 




n=l 



In (1 — x), 
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and, therefore, 



C(2) = Li 2 (l) 



In (1 — x) 



dx 



x 



dx dy 
1 — xy ' 



(4) 



o □ 
where □ = {(x,y) : < x < 1, < y < 1} is the unit square. Let us note 



dx dy 



+ 



1 — xy J J 1 + xy 

n 



dx dy 



dx dy 
1 — x 2 y 2 



and 



dx dy 



dxdy 1 f f dx dy 



1 — xy J J 1 + xy 2 J J 1 — xy' 

□ n 
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where the last equation follows from 



2xy 1 
— -dx dy = — 



1 — x 2 y 2 



dx 2 dy 2 



2JJ 1 — x 2 y 2 2JJ 1 — xy 
□ n 



dx dy 



It follows from equations ([5]) and (J6]) that 



c(2) ^//l 



□ 



dx dy 
— x 2 y 2 



(7) 



Now let us make the Beukers-Kolk-Calabi's magic change of variables in this 
two-dimensional integral [S] 



sin u sin v 
x = , y = 



COSf 



COSM 
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with the Jacobi determinant 



d(x,y) 
d(u, v) 



cos u sin v sin u 



sin u sin v cos v 



sin 2 u sin 2 v 

COS 2 f COS 2 u 



1 2 2 

1 — x y . 



Then miraculously 



C(2) = ^ J J dudv = ^Area(A), 



(9) 



where A is the image of the unit square □ under the transformation (x, y) — > 
(u,v). It is easy to realize that A is the isosceles right triangle A = {(u,v) : 
u>0,v>0,u + v<it/2} and, therefore, 

w . 4 1 /7T\ 2 7T 2 , , 

^ = 32(2) =T (10) 

"Beautiful - even more so, as the same method of proof extends to the 
computation of ((2k) in terms of a 2 k- dimensional integral, for all k > 1" 
[IT] . However, before considering the general case, we check whether the 
trick works for £(3). 



3 Evaluation of ((3) 

In the case of C(3), we begin with trilogarithm 

oo n 

^) = £^, (ii) 

n=l 

and using 

X 

d r ■ i \ t ■ i \ /* In (1 — J/) 7 
x— Li 3 (x) = Li 2 {x) = - / dy, 

dx J y 



we get 



i. X 

. dx f In (1 — y) , , . 

U3) - U,( 1 )= I — J -^—JLLdy. (12) 



o o 



But 



-I f^-y) dv = _[^-^ dz= f dz /_*_, 

y J a;^ J J I- xyz 

o o oo 

and finally 
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where D 3 = {(x, y,z) :0<x<l,0<y<l,0<^<l}is the unit cube. 
By the similar trick as before, we can transform f I X 3 p into the integral 

« 3 > = ?///i^ ("» 

n 3 

and here the analogy with the previous case ends, unfortunately, because the 
generalization of the Beukers-Kolk-Calabi change of variables does not lead 
in this case to the simple integral. However, it is interesting to note that the 
hyperbolic version of this change of variables 

sinh it sinhf sinht/; 

x = — ii > V = — r > z= — r (15) 
cosh v cosh w cosh u 

does indeed produce an interesting result 

C(3) = ^ JJJ dudvdw=*-Yo\{U 3 ) ) (16) 

where U 3 is a complicated 3-dimensional shape defined by the inequalities 

u > 0, v > 0, w > 0, sinhw < coshw, sinhf < cosh if, sinhw < coshw. 

Unfortunately, unlike the previous case, there is no obvious simple way to 
calculate the volume of U 3 . 

However, there is a second way to convert the integral (Tl2l) for ((3) in 
which the Beukers-Kolk-Calabi change of variables still plays a helpful role. 
We begin with the identity 

m = _jte hM^A iy= _ m<^v) dxdV: (17) 

J x J y JJxy 

D 

where the domain of the 2-dimensional integration is the triangle D = 
{(x,y) : x > 0,y > 0,y < x}. Interchanging the order of x and y inte- 
grations in the evaluation of the 2-dimensional integral (fT7|) . we get 

J y J x 

y 



5 



which can be transformed further as follows 

1 11 

/In ( 1 — yj\ny f f dx f f In y 
dy = - / In y dy / = - / / dx dy, 
y J J l-xy J J l-xy 

□ 

or in a more symmetrical form 

«3> = -5//r^*- ^ 



Note that 

2xy ln 2 ( f dxdy = \ lf^fWdtf = \ ff^dxdy. 
1 — x 2 y 2 4 J J 1 — x 2 y 2 4 J J 1 — xy 

□ □ □ 

Therefore, we can modify §J5§ and accordingly and using them transform 
(fig) into 

4 /Y ln(xy) 
7 77 



= T^adxdy. (19) 



At this point we can use Beukers-Kolk-Calabi change of variables (jSj) in ( TT9"j) 
and result we get 

C(3) = — - J J In (tan W tanw) dudv = — ^ In (tan it) dudv. (20) 

A A 

But this equation indicates that 

7I-/2 7T/2-U 7r/2 

C(3) = — - y dw In (tan u) J dv = — — J (^— — ixj In (tanit) du, 
o oo 

which after substitution x = | — w becomes 

tt/2 tt/2 



C(3) = — - / x In (cot x) dx — - / sin (tan a;) dx. (21) 



o 
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But 



tt/2 tt/2 

J In (tan x) dx = — J In (cot u) du = — J In (tan u) du = 0, 

tt/2 

which enables to rewrite (1211) as follows 

tt/2 tt/2 

C(3) 

o o 
and after integration by parts and rescaling x — » x/2 we end with 



tt/2 tt/2 

8 f / 7r \ . 8/",/ . d f x 2 tt\, 

- / x In tanx ax = - / In tanx — x ax, 

7 J V 4/ v 7 77 V ; dxV2 A ) ' 



C(3) = i/^& (22) 
7 ./ sin x 



o 



This is certainly an interesting result. Note that until quite recently very 
few definite integrals of this kind, involving cosecant or secant functions, 
were known and present in standard tables of integrals [121 D2J HI]- In fact 
( |22|) is a special case of the more general result [10] which we are going now 
to establish. 



4 The general case of £(2n) 

The evaluation of ((2) can be straightforwardly generalized. The polyloga- 
rithm function 

X) x" 



Li s {x) = Y,-s (23) 

obeys 



n=l 



x-^-Li s (x) = Lz s _i(x), 
dx 



and hence 



X 



Li s (x)= j^^-dy. (24) 

y 
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Repeated application of this identity allows to write 

1 XI X-n-2 

C(n) = Li n {l) = [*E±[*El... [ [- In (1 - x n ^)\ . (25) 



X\ J X 2 J Xn-1 



After rescaling 

x\ = x 2 = xiy 2 , x 3 = X 2 ?/3, • • • , x n -i = x n - 2 y n -i = y\Vi • • • Vn-l, 
and using 

i 

dy n 1 



o 



1 - Vm ■••Vn 2/l2/2 ■ ■ ■ Vn-l 



ln(l — 2/12/2 • • • 2/ 



• "Vn-l), 



we get 

C(n)= / ^■■■ dy \ (26) 



1 - 2/12/2 ■■•Vn 

□n 

where □„ is n-dimensional unit hypercube. The analogs of (jSJ) and ([6]) are 

/f dx\ ■ ■ ■ dx n f f dx\ ■ ■ ■ dx n f f dx\ ■ ■ • dx n 

J 1 - x l ■ ■ ■ x n J J 1 + Xi ■ ■ ■ x n J J 1 - x\ ■ ■ ■ x\ 

n n n n n n 

and 

dx\ ■ ■ ■ dx n f f dx\ ■ ■ ■ dx n 1 f f dx\ ■ ■ ■ dx n 



X "JC^ ' ,1 ,1 1 I Ob ]_ * * ■ X 2 / / 1 0u i * 0l> n 

C?i I— In I— Ijt. 

from which it follows that ( 1261) is equivalent to 

2 n -1 f f dxi--- dx n 

□„ 

If we now make a change of variables that generalizes (jSJ), namely 

sin Mi sinu 2 sinu n _i sinw n 
Xi = , x 2 = , . . . , x n _i = , x n = . (28) 

COS 1*2 COS COSM n COS Mi 
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we, in general, encounter a problem because the Jacobian of (125]) is jSJ E] 



d(xi,. ..,x n ) 
d(m, ...,u n ) 



1 (1) x 1 x 2 • • • x n , 



and, therefore, only for even n we will get a "simple" integral. For the 
hyperbolic version of (1281) . 

sinhf! sinht> 2 sinht> n _i sinht> n 

cosh v 2 coshu 3 coshf n cosh^i 

the Jacobian has a "right" form 

u{Xi, . . . , X n ) 22 2 

d(t>i, ...,v n ) 

and we get 



C(») = ■■■ d Vl ---dv n = — — Vol n (U n ). (30) 



However, the figure U n has a complicated shape and it is not altogether 
clear how to calculate its n-dimensional volume Vol n (£/„). Therefore, for 
a moment, we concentrate on the even values of n for which (|28l) works 
perfectly well and leads to 

C(2n) = ^— ^ J ■■■ J d Ul ---du n = Vol 2n (A 2 „), (31) 

where A n is a n-dimensional polytope defined through the inequalities 

A n = j(«i, . . . , u n ) : Ui > 0, + u i+ i < ^| . (32) 

It is assumed in (1321) that Ui are indexed cyclically (mod n) and therefore 

There exists an elegant method due to Elkies jH] how to calculate the 
n- volume of A n (earlier calculations of this type can be found in [15]). Ob- 
viously 

/7T\ n 

Vol„(A n ) = (-J Vol n (5 n ), (33) 



where Vol n (<5 n ) is the n-dimensional volume of the rescaled polytope 

S n = {(ui, ...,u n ):Ui> 0, Ui + m+i < 1} . (34) 

If we introduce the characteristic function Ki(u, v) of the isosceles right tri- 
angle {(u, v) : u, v > 0, u + v < 1} that is 1 inside the triangle and outside 
of it, then [3] 

1 1 n 11 

Vol n (5 n ) = J . . . J Y^Ki(ui,Ui + i)dui. . .du n = J du x J du 2 K 1 (ui,u 2 ) ■ ■ ■ 



o o 
i 







J du n -iKi(u n -2,u n -i) J du n Ki(u n -i,u n ) Ki(u n ,Ux). (35) 



Let us note that Ki(u,v) can be interpreted [9] as the kernel of the linear 
operator T on the Hilbert space L 2 (0, 1), defined as follows 

1 1— u 

(ff)(u) = J K 1 (u,v)f(v)dv= J f(v)dv. (36) 



Then (1351) shows that Vol n (5 n ) equals just to the trace of the operator T n : 

i 

Vo\ n (5 n ) = J K n (wi,Mi) dut, (37) 
o 

whose kernel K n (u,v) obeys the recurrence relation 

i 

m 



Surprisingly, we can find a simple enough solution of this recurrence relation. 



i 

K M = fK l{ u,u l )K^ ( u u v)iu l . 



Namely, 



KoJu,v) 



-l) n -„ -T X 



^2n-l ( — 7^ — ) + E 2n -i 



(2n - 1)! 
u — v 



E 2 n-i [ — t: — ) + E 2n -i 



2 

v — it 



9{u - v) + 
6{v - u) 



(39) 
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and 



2n 



2n 



K 2n +i(u,v 

1 — U + V 
2 

1 — U + V 



p2n-l 

-ir-^-T x 



+ E 2n 



-3 



2/i 



(2n)! 
1 — it — f 
2 

u + v — 1 



6(1 - u-v)+ 
9(u + v - 1 



(40) 



In these formulas E n (x) are the Euler polynomials [H] and is the Heav- 
iside step function 

r i, i/ x > o, 



h if x = 0, 



0, i/ x < 0. 

After they are guessed, it is quite straightforward to prove (J39l) and (1401 by 
induction using the recurrence relation (138j) and the following properties of 
the Euler polynomials 



^-E n {x) = nE^x), E n (l -x) = (-l) n E n (x). 
ax 



(41) 



In particular, after rather lengthy but straightforward integration we get 



i-tt 



K 2n +i(ui, v)dui = K 2n +2(u, v) - X, 



where 



X 



-1) 



n+l 



>2n 



But 



(2n + l)\ 

1 - v 



2n+l 



1 



1+V 

2 

1 + u 



2n+l 



1 - V 



2 2 

and the second identity of (T4T]) then implies that X = 0. 

Therefore the only relevant question is how f )39|) and f l40|) were guessed. 
Maybe the best way to explain the "method" used is to refer to the problem 
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13 from the aforementioned book Pp. To demonstrate the cardinal differ- 
ence between the ways problems are posed and solved by physicists and by 
mathematicians, Arnold provides the following problem for children: 

"On a bookshelf there are two volumes of Pushkin 's poetry. The thickness 
of the pages of each volume is 2 cm and that of each cover 2 mm. A worm 
holes through from the first page of the first volume to the last page of the 
second, along the normal direction to the pages. What distance did it cover?" 

Usually kids have no problems to find the unexpected correct answer, 4 
mm, in contrast to adults. For example, the editors of the highly respectable 
physics journal initially corrected the text of the problem itself into: "from 
the last page of first volume to the first page of the second" to "match" the 
answer given by Arnold [U [T7] . The secret of kids lies in the experimental 
method used by them: they simple go to the shelf and see how the first page 
of the first volume and the last page of the second are situated with respect 
to each other. 

The method that led to fl39|) and fj40l) was exactly of this kind: we simply 
calculated a number of explicit expressions for K n (u, v) using ( 13"8j) and tried 
to locate regularities in this expressions. 

Having ( 139|) at our disposal, it is easy to calculate the integral in (I37j) . 
Namely, because 

o2n-2 

K 2n (u,u) = (-1)" — [E 2n ^{u) + E 2n ^(0)] , (42) 

(2n — 1)1 



and 



we get 



/o2n-2 
K n (u,u)du=(-l) n - — ^2n-l(0), (44) 
(2n - 1)! 

o 

(note that E 2n (0) = E 2n {l) = 0.) But _E 2 n-i(0) can be expressed trough the 
Bernoulli numbers 

£ 2n _i(0) = -|-(2 2 "-l)B 2n , (45) 

and combining ( 13T1) . ( |33l) . ( )44l) and f j4"5]) . we finally reproduce the celebrated 
formula 

o2n-l 

a2n) = (-ir +1 —^B 2n . (46) 
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5 The general case of ((2n + 1) 

The evaluation of £(3) can be also generalized straightforwardly. We have 
C(n) = / Un - l{Xl) d Xl = f** J Ltn ~ 2{X2) dx 2 = [[ Un -^ d Xl dx 2 . 

J X X J Xi J X 2 J J Xi% 2 

D 

Interchanging the order of integrations in the two-dimensional integral, we 
get 

C(n)= j^^dx 2 [d_x 1 = _fln( X2 )L in . 2 (x 2 ) d ^ 
J x 2 J x x J x 2 

x 2 

Now we can repeatedly apply the recurrence relation ( 1241) . accompanied with 
Lix(x) = — In (1 — x) at the last step, and transform (j4"Tj) into 

C(n)= l^dxxl^... I ^± f ln ^~ X ^dx n . 2 , 

J Xx J X 2 J X n _ 3 J X n - 2 



which after rescaling 

x 2 = xxy 2 , x 3 = x 2 y 3 = Xl y 2 y 3 , x n _ 2 = x n _ 3 y n _ 2 = x x y 2 ■ ■ ■ y n -2, 
takes the form 

C(n)= f^dxx f ln{1 - Xm --- yn - 2) dy^ 2 . (48) 



Then the relation 
i 

dy n -i In (1 - xxy 2 ■ ■ ■ y n - 2 ) 



o 



1 - x x y 2 ■ ■ ■ y n -x ViV2---Vn-2 



shows that (j4"8"|) is equivalent to the (n — l)-dimensional integral 

/[ In 
••• / dxx ■ ■ -dxn-x- (49) 
J l-Xx---X n -l 
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As in the previous case, (j49p can be further transformed into 



C(") 



2 n - 1 



Inxi 



1 x x A n _l 



doc ~y d/CCyi — ]^ . 



or, in the more symmetrical way, 

2 n 1 
C(n) = - 



2"-lfi-l 



^ 1 2 — 1 dxi ■ ■ ■ dx n -i. (50) 



Let us now assume that n is odd and apply the Beukers-Kolk-Calabi change 
of variables (125)) to the integral (JSDJ). As the result, we get 

1 2 2n+1 /" f 
((2n + 1) = - — 22n+1 - ^ ■■■ In [tan («i) • • • tan (u 2n )] dui • • • du 2n , 



which is the same as 



C(2n+1) 



)2n+l 



2 2n+l _ J 



In [tan (tti)] diti ■ ■ - du 



In- 



By rescaling variables, we can go from the polytope A 2n to the polytope 5 2n 
in this 2n-dimensional integral and get 



C(2n+r 



22n+l 



22n+l _ 1 V 2 



7T 



2n 



In 



7T 

tan ( ui — 



dtti ■ ■ ■ dU2n- (51) 



Using the kernel K2 n (u, v), we can reduce the evaluation of ( 15"TT) to the eval- 
uation of the following one-dimensional integral 



C(2n + 1) 



2vr 



2/i 



2 2n+l _ ]_ 



But 



In 



7T 



In 









tan ( 













7T y 




cot ( 


—u 






2 /J 





K 2n (u,u) du. (52) 



7T 

tan ( —V 
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which enables to rewrite (j52p as 
C(2n+1) = 



7T 



2» 



j In 




tan ( K ^ u ^j 







2 2n+i _ i I — l — ' v 2 j ^ 2 "^' M ) ~ ~ u ' 1 - M )] 



(53) 

However, from (jl2j) and (l4"3"j) we have (recall that £? 2 n-i(l — w) = — -E , 2 n-i(^)) 

K 2 n(M, - ^2n(l - U, 1 - «) = 



o2n-l J 



(2ra)! 

and the straightforward integration by parts in ([SSD yields finally the result 

i 



C(2n + 1) 



;-i)»tt 



n „2n+l 



4[1 _2-(2«+i)] (2 n )l 7 sin(7rn 



EoJu) 



■ du. 



(54) 



This is exactly the integral representation for ((2n + 1) found in [10]. Our 
earlier result ( l22p for ((3) is just a special case of this more general formula. 



6 concluding remarks: £(2) and amoebas 

It remains to clarify the origin of the Beukers-Kolk-Calabi's highly non-trivial 
miraculous change of variables f )28|) . Maybe an interesting observation due 
to Passare [7] that £(2) is related to the amoeba of the polynomial I — Z1 — Z2 
gives a clue. 

Amoebas are fascinating objects in complex geometry (THJ [19] . They are 
defined as follows [20]. For a Laurent polynomial P(zi, . . . , z n ), let Zp denote 
the zero locus of P(zi, . . . , z n ) in (C\{0}) n defined by P(z±, . . . , z n ) = 0. The 
amoeba A(P) of the Laurent polynomial P(z\, . . . ,z n ) is the image of the 
complex hypersurface Zp under the map 

Log: (C\{0}) n M n 

defined through 

(z 1 ,...,z n ) -> (ln|zi|,...,ln|z n |). 
Let us find the amoeba of the following Laurent polynomial 

P(zi, z 2 ) = Zi - zl x -i(z 2 - z^ 1 ) . (55) 
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Zl = e u e*+", z 2 = e v e- i ^, 



Taking 

we find that the zero locus of the polynomial fl55|) is determined by conditions 
cos <p u sinhw = sint^ coshf, sin0 u coshw = cos0„ sinh v. 

If we rewrite these conditions as follows 

sinhf sin0 u sinhw sin0„ 

x = — TT = 1 - ' y = — r — — > 56 

cosnu cos^ coshf cos0 u 

we immediately recognize the Beukers-Kolk-Calabi substitution ([8]) and its 
hyperbolic version with the only difference that in (jHJ) we had < x,y < 1. 
However, from f )56l) we get 

cos 2 <p u = — ^ , cos 2 4> v = — — |— , (57) 
1 — x z y z 1 — x z y z 

and 

2 1 + ?/ 2 , 2 1+a: 2 

cosh m = — , cosh u = — -. (58) 

1 — x z y z 1 — ary 2 

It is clear from fl57|) and (158]) that we should have 

x 2 < 1, y 2 < 1. 

Therefore, the amoeba A(P) is given by relations 

* i \ f , . sinh it sinh t> 1 . . 

A(P) = \ (u,v) : -1 < — — < 1, -1 < — — < 1 , 59 
[ coshf coshw J 

and the hyperbolic version of the Beukers-Kolk-Calabi change of variables 
(JHj) transforms the unite square □ into the one-quarter of the amoeba (|59|) . 
Then the analog of ([9]) indicates that £(2) equals to the one-third of the area 
of this amoeba. 

As we see, the hyperbolic version of the Beukers-Kolk-Calabi change of 
variables seems more fundamental and arises quite naturally in the context of 
the amoeba fl59|) . Trigonometric version of it then is just an area-preserving 
transition from the "radial" coordinates (u, v) to the "angular" ones (<p u , <f> v ). 

One more amoeba related to £(2) was found in [7]. Although the corre- 
sponding amoeba .4.(1 — z\ — z 2 ) looks different from the amoeba ( T59l) . they 
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do have the same area. The trigonometric change of variables used by Pas- 
sare in [7] is also different from (jSJ) but also leads to simple calculation of the 
area of A(l — Z\ — z-i) and hence C(2)- Of course it will be very interesting 
to generalize this mysterious relations between ((n) and amoebas for n > 2 
and finally disentangle the mystery. I'm afraid, however, that this game is 
already not for kids under fifteen. 
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